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Q\ . Abstract 

The electromagnetic field inside an isotropic optically active medium is quantized on the 
base of Fedorov's model for optical activity. The modified photon propagator is derived. Using 
I this result it is shown that the QED correction to the electron mass contains chiral terms which 

r ^ . nevertheless are parity invariant. The spontaneous emission rate of a two-level atom is shown 

to increase in the medium. The spontaneously emitted radiation is partially polarized. 

t> ' 

PACS: 12.20.-m,33.55.Ad 



Introduction 



- 1—^ 

X 



In recent years there has been a growing attention to quantum electrodynamics (QED) in dielectric media 
jl], ||, This field is of interest for several reasons. One motivation is the study of squeezed states 

' of light connected with the dielectric ||. Another one is the more accurate description of Casimir forces 

between conducting plates which are special cases of dielectrics (see, e.g., Ref. Q and references therein). 
Despite these efforts the case of optically active media, for which (in most models) the electric displacement 
vector depends also on derivatives of the electric field, was not often studied. To the author's knowledge there 
is only a work of Eimerl in which the photon self energy is used to give predictions for the parameters 
Oh, describing the medium, and a paper of Woolley || where a nonrelativistic QED calculation for the molecular 

1 magnitudes leading to optical activity is presented. 

The influence of the quantized electric field on quantum matter in the presence of an optically active 
medium seems not to be discussed in the literature. In this paper it is considered how mass renormalization 
for the electron and spontaneous emission will be modified by the medium. Compared to ordinary dielectrics 
new effects mainly arise because of the different behaviour of left and right circularly polarized field modes. 
Since these polarization modes are the eigenmodes of the spin operator many new effects depending on the 
spin and the angular momentum of atoms or electrons can be imagined. This is the case for the self energy of 
the electron, for instance, which will be shown to include chiral terms. The spontaneous emission of photons, 
which in the approximation made here does not include spin or angular momentum of the atom, leads to 
partially polarized light as will be demonstrated later. Beside the examples studied here the Casimir effect 
for optically active media is under investigation. It should be mentioned, however, that the effects produced 
by the medium are in general small since the change in the refractive index in optically active media is small 
compared to the change in ordinary dielectrics. 

Quantization of the electromagnetic field 

The starting point of the calculations are the Maxwell equations 

V-D = p,\7xH-D = j 
\7-B = 0,VxE + B = 

together with Fedorov's model for optical activity B 



D = e(E + /3V x E) 

B = n{H + /3V x H) (2) 
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where only the case of an isotropic medium is considered. For the sake of simplicity e and /x, which are not 
related to the optical activity, will be replaced in this paper by their vacuum value, e = eq, [i — /io- Fedorov's 
model was chosen to preserve gauge invariance (for a compilation of other models see Ref. ]h|). The reason 
for this is that the nontrivial part of the material equations (^) depends on E and B only through their curl. 
This implies that only the transverse electromagnetic field, which is not affected by gauge transformations, 
will be modified by the medium. 

Switching to the gauge potentials (f>, A and imposing the Lorentz gauge condition on them one easily 
deduces the field equations 

□0 = 

uA = ~4^ 2 { 2V xl + /JVxVxI} (3) 



in absence of free charges p or sources j. The d'Alembertian □ is defined by d 2 /c 2 — A. To quantize these 
fields the method of normal mode decomposition will be used. To do so one first shows that the scalar 
products 



(</>>) := / d 3 x{<j>^-<f>^} 



(A, A) := J d 3 x^A ■ A - A^ ■ A + 2j8[(V x A) ■ A - $ ■ (V x 1')] + 

B 2 [(V x A) • (V x A) - (V x ft) • (V x 1')] } (4) 

are conserved for solutions <f>, ip and A, A' of Eqs. (|^). The general solution of the field equations can be 
found by Fourier transformation. It is given by 



(2ir) 3 e cJ \ V2k . /^(l - /3fc) J2k£.(l+Pk) 



Here efj are the complex polarization vectors defined by e*3 := k/k and ie^ x e± = ±e±. Although the 
calculation is not Lorentz invariant it is convenient to adopt the covariant notation a ■ b := a^b^ = a°b° — a-b 
and x° := ct. Throughout the paper the summation convention over repeated indices is adopted. Greek 
indices are running from to 3 and latin ones from 1 to 3. The various four-wavevectors are defined by 
fc' 1 := (fc, k) and fc^ := , k) where 

4 := k ■ (6) 

y/l±0k- (3 2 k 2 

Here one can see the main effect of the presence of the optically active medium. The dispersion relation 
between frequency and wavevector of the left and right circularly polarized modes (±) has been changed from 
k° = k to fc° = k±. Since k± diverges for plus-modes with (3k = 1 + \/2 and minus-modes with (3k = y2 — 1 
it is clear that the theory describes physics only for electromagnetic waves with (3k <C 1. This should not 
be regarded as a shortcoming since in real media optical activity appears only well below the high energy 
physics regime. The introduction of a cutoff in momentum space is therefore physically meaningful. 

It should be mentioned that the dispersion relation k° = k± is a good description only for a certain 
part of the low energy frequency range. This is a consequence of describing optical activity within a local 
field theory. Despite this fact local theories can still be a useful tool estimate the order of magnitude of the 
medium's influence on quantum matter. 

In Eq. (||) it also becomes apparent that only the transverse field modes (±) feel the presence of the 
medium in Fedorov's model. Hence no gauge transformation can alter the effect of the medium on the 
Maxwell field. 
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The quantization of the Maxwell field is now an easy task. One demands the usual Lorentz gauge 
commutation relations 

[a a (k),al{k')] = 5(k-k , )5 af3 (7) 

with a, (3 = 0, +, — , 3 and — <5oo = $33 = <5++ = 8 = 1 as well as 8 a p — for any other combination of a 

and (3. In addition one has to impose the Gupta-Bleuler condition d^A^+'lip) = on the physical states 
\t/>) (see, e.g., Ref. 0). 

The self energy of the electron 

To calculate the first order mass correction for the electron the Feynman propagator 

D^(x -y) = -i6{x° - y Q )[A^ + \x),A^-\y)] - i9{y° - x°)[A"^ (y), v4 m( ~) (x)] (8) 

of the photon field is needed. From Eq. (||) it is clear that only the spatial part of the propagator is affected 
by the medium. It is useful to rewrite D 1J into the form D l J ac + D l _l + D l _l where D l J ac stands for the vacuum 
contribution and D± describes the changes of the propagation of the left and right circularly polarized 
modes induced by the medium. Using standard methods (including an integration over k° with the residue 
formula and the use of e+(— ( — k) — e+(fc)e+(&) and a similar equation for e_) the expression for D l ± 
in momentum space is found to be 

D ± (kn = ~/ ±{k) ^ {k) { (1 T /?fc) 2 [(fc°) 2 - (4 ) 2 + ie] ~ (fc°) 2 -fc 2 + *e} ' (9) 
Having found the photon propagator one is ready to calculate the the first order on-shell mass correction 

8m = J d 4 kD^(k) la S F (g - khp (10) 

with Sf(q) — [l^lfi ~ A* + * e ] 1 being the Fermion Propagator, hq^ is the external momentum of the 
electron and \i := mc/h is its Compton wavevector. If j^qn stands on the right of dm it can be replaced 
by /i since the external electron spinor ip fulfills the free Dirac equation (7^9^ — n)ip — 0. Again it is 
convenient to write 8m in the form Sm vac + Sm + + 8m^ where every term corresponds to the respective 
part of the photon propagator. Then 5m± can be calculated in a standard manner by making a Feynman 
parametrization and performing the k° integral with the residue formula. Using a^j^jj = — 7 :) a A1 7 A1 + 2a j 
and —(i/2k)eiji'f l j : >k l k a '-f a — — k°j5jik l /k + £7075 as well as 

elel* = i{^--ifc < FT^^} (11) 



(e 123 = — 1) one arrives at 

=2 r 1 



8m± ■ 



l ^/^{^ (1 T ^3} 



(27r) 3 4e c 2 

{ ± 7 o7 5 fc + Tzq°l57ij ~ (<f + fc*bi - zq° lo + ^ 7 ,(fc • q) T ie^'Sfc^} (12) 



with w± := J z 2 /j, 2 + (k — zq) 2 + (1 — z)((k±) 2 — k 2 ) and where w is defined by the same expression if k± is 
set equal to k. After fixing the direction of q to be parallel to the 3-axis the angular part of the integration 
can be performed as usual. In the remaining expression a high-frequency cutoff k is introduced for the dk 
integral. As explained above this cutoff is needed since optical activity is only present in the low-energy 
regime. Since k has to be chosen so that (3k <C 1 is fulfilled it is sufficient to calculate 8m only to first order 
in (3. In this case the remaining integrations over z and k lead to closed expressions for 8m which contain 
various polynomials and logarithms. For all allowed external momenta (q(3 < 1) the lowest order Taylor 
expansion 

e 2 (3k 3 



Sm_ = 



9ir 2 Soc 2 [i 2 

gives an excellent approximation to these functions 



{ - 3^7075 + 97573} + 0{(3 2 ) + 0(q 3 ) + 0{k a ) (13) 
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This result deserves some comments. First, it is obvious that 5m does not only depend on q^q^ as in 
the vacuum but also on the three- vector q. This implies that for on-shell amplitudes with q^q^ = /i 2 the 
propagation depends on the momentum of the electron in a non-trivial way. The reason is that the rest 
frame of the medium introduces a new time-like vector £ M so that 5m can depend on both q^q^ and q^C,^- 

The second comment to Eq. ( |l3| ) concerns the fact that 5m+ + 5m^. depends on the spin of the electron 
only through the chiral operators 7075 and 757^ which are known to be a pseudo-scalar and pseudo- vector, 
respectively. Nevertheless, 5m behaves as a scalar under a parity transformation if one takes into account 
the transformation properties of the medium. For Fedorov's model the parameter [3 is a pseudo-scalar |j, [l^] 
so that 5m indeed has the right transformation properties. 

The physical origin of the chiral mass corrections is hidden in the properties of the medium. As explained 
in Refs. jl3], the self energy of a particle is a consequence of the radiation reaction and not of the vacuum 
fluctuations of the electromagnetic field. Hence it is produced by the held emanated from the electron that 
interacts with the medium and turns back to the electron. The chiral nature of the self energy is caused by 
the different interaction of the medium with left and right circularly polarized light. This implies that held 
modes with different total angular momentum propagate differently (as can be seen in Eq. ([)[)) and the back 
reaction on the electron depends on its total angular momentum. 

To estimate the magnitude of 5m + + (5m_ a specific value for (3 is needed. This can be done by applying 
the method of Ref. |l(J to Fedorov's model. The result is that (3 is given by X(n + — rt_)/(47r) where A is 
the wavelegth of the light beam and n± are the refractive indices for left and right circularly polarized light. 
Adopting the values for quartz (that, however, is not an isotropic optically active medium) at a wavelength 
of 762 nm one can infer that /3 is given by 3.6 x 10~ 12 m. Then a very large value for k is, e.g., k = 10 11 
m _1 . Inserting these values into Eq. ( |l3| ) gives (<5m+ + 5mJ)jm i=a 10~ 6 , the effect of the medium on the 
electron mass is very small. 

Spontaneous emission of a two-level atom 

Spontaneous emission of a photon during an atomic transition and the exponential decay of the excitation 
probability arise from the destructive interference between the transition amplitudes to all photon states. 
The calculation of the decay factor for a two- level atom with excited state |e) and ground state \g) is a 
standard task in quantum optics and can be found in many textbooks (see, eg. Ref. ||). In the Markov 
approach the Hcisenberg equations of motion of the atomic density matrix and the annihilation operators of 
the radiation field is solved. In this paper the Hamiltonian 

H = H A + H e . m . - d- E x (\e)(g\ + \g)(e\) (14) 

is used, where Ha '.= -E e |e)(e| + E g \g)(g\ contains the internal energy of the atom and H e . m describes the 
time evolution of the free radiation field, d :— (e\qx\g) is the dipole moment of the atomic transition and 
E 1 - = —dtA ± is the transversal part of the electric field. In Fedorov's model E 1 - can be derived from Eq. 
(|J) and the free evolution of the Maxwell field is governed by the Hamiltonian 

H e . m . = f d?k hck° x a{(k)a x (k) (15) 

where the ground state energy was removed. Since the atom couples only to the transversal part of the 
radiation field the longitudinal and scalar part was neglected in H e . m .. 

The essential step of the Markov approach is the Markovian approximation in the formal solution for the 
annihilation operators (see, e.g., Ref. ||). Following these steps the decay rate of the excited state is found 
to be 

~f = 2ir[d 3 kJ2\Cx{k)\ 2 6(ck Q x ~uu ). (16) 

Here wo '■— (E e — E g )/Ti is the atomic transition frequency and C\{k) are the expansion coefficients of the 
electric field operator defined by d ■ E =: Ti J d 3 kJ2x=±(^>~(^) a x(^) + h.c). The 5 distribution guarantees 
energy conservation. It also implies that the Fedorov model as well as any other local theory of optical activity 
can predict the spontaneous emission rate much better than the mass correction since only field modes with 
a certain frequency can contribute to it. The variation of the parameter (3 with the field frequency, which can 
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only be included in nonlocal theories, does therefore not alter the results. After evaluation of the integrals 
in Eq. (fLq) the contribution j± of the left and right circularly polarized modes turn out to be 



7o 4 



(17) 



7± = 



2 k 5 (lT[3k ± ) 2 (l±Pk±) 



with fco := uio/c. The wavevector belonging to modes with energy Twjq is given by k± :— &o(\/l + 2/3 2 fcp ± 
(3k )/(l + f3 2 k^). The decay rate in free space is denoted by 70 := (Pk^/ {Zi:Tieo). To a very good approxi- 
mation the total decay rate 7 = 7+ + 7- is then given by 



This result shows that an isotropic optically active medium always increases the spontaneous emission rate. 
This is due to the fact that the decay factor, being the inverse of the lifetime, has to transform as a scalar 
under parity transformations. Hence only even powers of /? can occur in 7. That the coefficient of 1 is 
positive can be understood by looking at j± as a function of k±. Since k± = ko(l ± fiko) + 0((3 3 kg) we 
need only to consider the lowest order term in k± depending on (3. Then k + and fc_ are shifted by the same 
amount from their vacuum value, but with opposite sign. As is well known Eq. (Hq) states that most of the 
/c-dependence of 7 is caused by the density of the field modes (see, e.g., Ref. ||). This density increases 
more if k is shifted to a higher value (k+) than it decreases if k is shifted by the same amount to a lower 
value (fc-). Thus the grow of 7+ is larger than the decrease of 7_ so that the total decay rate increases. If 
(3 is less than zero 7_ is larger than 7+. This argument also shows that the spontaneously emitted radiation 
is partially polarized since for (3 > 0, say, more left circularly than right circularly polarized photons are 
produced. 
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